NASA/TM- 207522 


y/y-^Y- T'^ 



_ AIAA 80-1383R 
Nonreflecting Far-Field Boundary Conditions for 
Unsteady Transonic Flow Computation 

D. Kwak 


Fteprinlad Irom 

HUM 


Volume 19. Number 1 1 . November 1981 , Page 1401 

This paper is declared a work o1 the U S. Government and therefore is in the public 
domain 


AMERICAN INSTITUTE OF AERONAUTICS AND ASTRONAUTICS • 1290 AVENUE OF THE AMERICAS • NEW YORK. NEW YORK. 10104 


VOL. 19, NO. 11, NOVEMBER 1981 


AIAA JOURNAL 


1401 


AIAA80-1393R 


Nonreflecting Far-Field Boundary Conditions 
for Unsteady Transonic Flow Computation 

D. Kwak* 

NASA Ames Research Center, Moffett Field, Calif. 


The approximate nonreflecting far>field boundary condition, as proposed by Engquist and Majda, is Im- 
plemented In the computer code LTRAN2. This code solves the Implicit finite-difference representation of the 
small-disturbance equations for unsteady transonic flows about airfoils. The nonreflecting boundary condition 
and the description of the algorithm for implementing these conditions in LTRAN2 are discussed. Various cases 
are computed and compared with results from the older, more conventional procedures. One concludes that the 
nonreflecting far-field boundary approximation allows the far-field boundary to be located closer to the airfoil; 
this permits a decrease in the computer time required to obtain the solution through the use of fewer mesh 
points. 


I. Introduction 

T O compute time-dependent flows over an infinite region 
using finite-difference procedures requires either that the 
problem be reduced to the one in a finite domain via a 
coordinate transformation, or that boundaries be placed at a 
finite distance from the body and the boundary conditions 
modified along them. When the solution at infinity is simple, 
a coordinate transformation technique can be used. However, 
for many practical problems, especially when the flow is 
singular at infinity, the coordinate transformation method is 
not applicable. In such cases, computational simulations of 
unsteady flows in an unbounded region are performed on 
grids with finite dimensions. The artificial-wall effect created 
by these grid far-field boundaries must therefore be 
minimized so as not to degrade the resulting numerical 
solution. 

Various approaches have been developed to reduce this 
spurious influence of far-field boundaries on interior 
solutions, some of which are listed below: 

1) It is common' to place these boundaries far enough 
(typically hundreds of chord lengths) from the airfoil such 
that reflected waves have been greatly reduced in amplitudes 
by numerical dissipation through the travel from these 
boundaries. Such a procedure naturally requires a large 
number of grid points. 

2) In a series of papers, Engquist and Majda developed a 
nonreflecting boundary condition procedure. By absorbing 
waves incident on the boundary in certain directions, they 
designed approximate local conditions for simulating the far- 
field effects along those boundaries. Their lowest order 
approximation can be interpreted as a simplified form of the 
characteristic equation. 

3) Bayliss and TurkeP derived a far-field boundary 
condition procedure based on an asymptotic expansion to any 
order in the reciprocal of the distance from the origin. The 
dependence of the far-field distance is shown in the work of 
Fung,^ where the linearized unsteady small-disturbance 
transonic flow equation is solved to get a far-field condition. 
This condition is equivalent to a vortex solution with a time 
lag for waves originating at the airfoil and propagating to the 
far-field boundary. 

4) A viscous damping method which absorbs incident 
waves on the far-field boundaries was tried by Bushby and 
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Timpson."^ This method, however, requires a considerable 
number of grid points near the far-field boundary. 

5) Guderley* (see also Krupp and Cole’) has shown that in 
a limited number of cases, far-field conditions can be derived 
analytically to investigate basic mathematical concepts. 
However, an extension of this to more general nonlinear 
problems is not clear at this point. 

Unsteady transonic flow regions are encountered in the 
flight of many existing flight vehicles. Typical examples 
include a fluttering airfoil and spinning helicopter rotors. 
For unsteady small-disturbance transonic flows, Ballhaus and 
Goorjian' have developed a time-accurate implicit finite- 
difference computer code, LTRAN2. Their boundary con- 
ditions at the lop, bottom, and upstream boundaries (see Fig. 
I) are perfectly reflecting conditions and arc the proper ones 
for boundaries infinitely far away. In the present work, 
approximate nonreflecting boundary conditions are employed 
in this code to bring the far-field boundaries closer to the 
airfoil. As illustrated in Fig. 1, waves originating from the 
airfoil propagate to the far-field boundaries. The old 
boundary which is placed at a large distance from the airfoil 
leaves the near-field solution unaffected by reflected waves 
from the outer boundary. With nonreflecting conditions, the 
new boundary can be located closer to the airfoil, thus 
reducing the area of computation and thereby increasing the 
degree of resolution within the computational domain. 

Among the approaches listed above, the approximate local 
condition procedure designed by Engquist and Majda was 
selected to be used in the present study. Their conditions are 
simple to implement and require only a modest change in the 
existing computer code; however, these conditions result in a 
marked improvement in the computational efficiency of that 
code. 

In Sec. II of this paper, the governing equation and the 
boundary condition procedure by Engquist and Majda are 
presented. In Sec. Ill, the finite-difference scheme for im- 
plementing these conditions in LTRAN2 is described. In Sec. 
IV. computed results using the old reflecting far-field 
boundary conditions and those using the new nonreflecting 
conditions are compared with large grid solutions. 

II. Governing Equation and Far-Field 
Boundary Conditions 

A. Governing Equation 

The unsteady, transonic small-disturbance equation for low 
reduced frequencies can be written as 

= + (la) 
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Fig. I Schematic of computalional domain: old reflecting boundary 
vs new nonreflecting boundary. 


where 

+ (lb) 


Here, 0 is the disturbance velocity potential, 6 the airfoil 
ihickness-lo-chord ratio, and the freestream Mach 
number. The reduced frequency is defined as k^ioc/U^ for 
an airfoil of chord length c executing some unsteady 
oscillatory motion of frequency w. The choice of the exponent 
m is somewhat arbitrary and is made to extend the Mach 
number range of the small-disturbance theory (a rather 
complete review is given in Ref. 10). The quantities x, y, /, and 
0 have been scaled by c, and respec- 

tively. In deriving Eq. (la), it is assumed that 

( 2 ) 


B. Far-Field Boundary Ct>ndltions 
In LTRAN2, Eq. (la) is solved for the flow about airfoils 
by an alternating-direction implicit (ADI) finite-difference 
algorithm. The existing far-field boundary conditions Im- 
posed in LTRAN2 are 


upstream: 

<t> = 0 


downstream: 


(3) 

top and bottom: 

<t>,=0 



These conditions are perfectly reflecting. In other words, all 
of the waves originating from the airfoil for unsteady 
problems are reflected back into the computational domain 
from the boundaries. Hence, the outer boundaries must be 
placed far enough from the airfoil in the original version of 
LTRAN2. 

Following Engquist and Majda,^ conditions are introduced 
to absorb at least a portion of the waves incident on the outer 
boundaries. For the upstream condition, B in Eq. (lb) i$ 
assumed to be locally constant; therefore, the analysis is 
performed on a linear equation. This approximation assumes 
that the flowfield is governed by a linear equation in the 
vicinity of a far-field boundary. By considering waves 
traveling left from the interior to the upstream boundary for 
subsonic freestream, a perfectly nonreflecting condition is 
derived. Since this condition requires information from the 
mathematical domain of dependence of that boundary, 


approximations are made to get the following local conditions 


upstream 


1st approx: 0, — 0=0 
A 

^ ^ A 1 

2nd approx: 0^, = (j 


(4) 

(5) 


The first approximate condition absorbs waves normally 
incident on the boundary, while higher approximations 
absorb portions of obliquely incident waves. 

Since waves travel infinitely fast in the downstream 
direction, the following consistency conditions seem ap- 
propriate (see Krupp and Cole’) 

downstream Tor 

= 0 (7) 

Considering waves traveling to the top and bottom 
boundaries, a nonreflecting condition is developed. Sub- 
sequently, local approximations to it are made, resulting in 

top and r 1st approx: 0^±VT^0j^ = O (8) 

bottom I 2nd approx: 0,^, +r^0^j, ±r/0^, =0 (9) 

Here, ± corresponds to top and bottom boundaries, and 

govern the absorption of waves obliquely incident to the 
boundaries. The values T;, are discussed in more detail by 
Engquist and Majda.^ 

Alternatively, a nonreflecting far-field condition can be 
designed by considering characteristic equations. Assuming 
5 = const in Eq. (lb) at a large distance from the airfoil, the 
characteristic equation can be written as 

+ = 0 (10) 
This equation is satisfied by the characteristic plane 

Hx,y,t) =r- {x/\TB^‘2^t/A) =const (I la) 


where 


r= {x/yfB,y) 
and 

0=/(f) (lib) 

is a solution to Eq. (la) which represents a traveling plane 
wave. From these, the following equation is derived by 
forcing linear combinations of the derivatives of 0 to be zero. 



A X 
B~r^ 



(12a) 


This equation can be used as a nonreflecting far-field 
boundary condition. Since this equation requires the 
evaluation of x/r and y/r along the boundary at each time 
step, a further simplification is made. In the upstream region, 
x^ - 00 with finite, so that Eq. (12a) becomes 


0,-(fi//1)0^ = (? (12b) 

For the top and bottom boundaries, y^ ±oo with Ijrl finite, 
so that Eq. (12a) becomes 

0,±/i(/I/\^)0,=O 

Replacing 0^ by 0^ using Eq. (lib), this equation can be 
written equivalently as 

0 ,, ± v^0 ,. = 0 


(12c) 
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These are the first approximate relations obtained by 
Engquist and Majda. Therefore, we may interpret Eqs. (4) 
and (8) as special forms of a far-field characteristic relation 
designed for simple, yet very efficient, applications. 


Ill . Finite-Di fference Scheme 

In LTRAN2, ' Eq. (la) is solved by a two-step procedure to 
advance from time step n to n -t- 1 


jf sweep: /I (At) 
y sweep: /l(AO 


The operators in Eqs. (13) and (14) are defined by the 
following equations 

(first order) 

(second order) 

fjj- ym, ^ 

= for + 


Here, Ar is the lime step and 7, t are the grid point indices in 
the X and directions. 

The first approximate nonreflecting far-field boundary 
conditions are readily incorporated into the above two- 
dimensional sweep procedure. For the x sweep, the upstream 
boundary condition [Eq. (12b)] becomes 

(A/) - ' («"/„., - 4>;. „ 2 .,) - - ' ( Ax) - ' 

j=l, f= /,..., ftnax 


Here 

B,y „2.t = - Mi ) - 1(7 + / ) (4>;. I.t - Ki'> 

x(x^+,-Xy) -'J. J-> 

The downstream boundary condition is 


^ymax.f - U ^jman - U ^ 


(16a) 


or equivalently, since from the previous 

iteration 


f= /..••.ftnax (16b) 

Here, ymax and ftnax arc maximum indices in x and y 
directions and represent downstream and top boundaries, 
respectively. 


Similarly, for the y sweep, the top and bottom boundary 
conditions become 

-y<- / ) =*= 

X (</>", t'//; - <!>"-!}- I,2)nxj-x,.,)=0 

y=2,...,ymax (IT) 

where 

^nyll2 +^'<) 

When j = 2, Eq. (17) requires aly = I . Therefore, fory = 2 we 
use Eq. (12b) also. Here, ihe upstream condition Eq. (15) is 
changed to 

(A/) (^;.v,zf-<^;.//y.f) - y^B^.u2,fA - ' (Ax) 

X I ( M /.f - <^;.r ) 1 = 

y=/, f=/, or fmax (*8) 

Then Eqs. (17) and (18) are solved simultaneously. 

The higher-order approximations can be applied in a 
similar way, however, at the expense of increased computing 
time. 


IV. Computed Results for Various Test Cases 

In the production version of LTRAN2, the default grid 
boundaries are located 857 chords from the airfoil in y 
direction and 200 chords from the leading edge, with 
smoothly stretched grid spacings in both directions (113, 97 
mesh in x,y directions). The large grid solutions in Figs. 2-7 
are obtained using this grid. However, depending on the 
particular problem, equivalent results can be obtained with 
smaller grids, and the minimum number of mesh points 
required is shown in the figures. Therefore, identical solutions 
are obtained when far-field boundaries are placed farther 
away than the minimum required. 

In computing the large grid solutions, two different sets of 
far-field boundary conditions are imposed independently, 
namely: 1) perfectly reflecting conditions as given by Eq. (3); 
and 2) the first approximate conditions as given by Eqs. (4), 
(6), and (8). These boundary conditions yield essentially 
identical results, i.e., within plottable accuracy. Since the 
large grid results can be duplicated using different com- 
binations of grids and boundary conditions, it is reasonable to 
assume that these are solutions of Eq. (I) with far-field 
boundaries at infinity. 

A. Step Change In Angle of Attack 

An impulsively started airfoil in plunging motion creates a 
pressure pulse. This is simulated numerically by a step change 
in the angle of attack. Due to the low-frequency nature of the 
governing equation, the lift and moment coefficients 
gradually increase after the initial change in the angle of 
attack. A problem such as this provides a good test case for 
investigating the inHuence of the computational far-field 
distances on the flow solution near the airfoil. 

Starting from the standard large grid, the far-field 
boundary distances are reduced by removing outer grid lines 
in the x (both upstream and downstream equally) and y 
directions. Both the old reflecting boundary conditions [Eq. 
(3)] and the new nonreflecting ones are tested with these 
reduced grids. In Figs. 2 and 3, the lift and moment coef- 
ficients for an NACA 64A006 airfoil are plotted on a lime 
scale in units of chord lengths of airfoil traveled. The indicial 
responses to a unit change in the angle of attack a are 
and as described in Ref. 11. and can be regarded as 
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b) MIDCHORD PITCHING MOMENT COEFFICIENT 

Fig. 2 Effect of far-field boundary distance on C(^ and vs / for 
step change in angle of attack. 


Cf/a and C,„/a in this case. By using the old conditions, the 
computed results start to deviate from the large grid solutions 
when the mesh gets smaller than 103x81 at =0.80 and 
0.85. However, applying the first approximate nonreflecting 
conditions, the mesh can be reduced to 93x71 at =0.80 
without significant deviations from the large grid results. 
Similar results are obtained at =0.85, as shown in Fig. 2. 

Since the x boundary can be brought to 3.8 chords from the 
leading edge using the first approximate condition, the second 
approximation was tried to bring the y boundary closer. 
However, there was practically no improvement to justify the 
additional computation required to implement the second 
approximation. In Table 1, the minimum mesh and the 
computational efficiency required to use the old and the first 
approximate boundary conditions are shown. Using the 
nonreflecting boundary conditions, the computational ef- 
ficiencies gained are 19 and at Afo,=0.80 and 0.85, 
respectively. 

B. Oscillating Flap Case 

Another test of the nonreflecting boundary conditions was 
performed for the case of an oscillating trailing-edge flap. The 
configuration consisted of an NACA 64A006 airfoil with a 
sinusoidally oscillating flap. Of the three types of possible 
shock-wave motions, type A (sinusoidal shock-wave motion) 
and type B (interrupted shock-wave motion) are shown in 
Figs. 4 and 5. In both cases, large grid results can be obtained 
with a grid as small as 93x77. Again the minimum x 
boundary distance is very small. As expected, the influence of 
the y boundary is very pronounced. By the time the y 
boundary distance is reduced to 7.6 chords from the airfoil, 
the old reflecting conditions cause significant deviation of the 
lift and the moment coefficient from the large grid solution. 


ASYMPTOTIC VALUES 



a) LIFT COEFFICIENT 


ASYMPTOTIC VALUES 



0 10 20 30 40 50 60 

TIME, CHORD LENGTHS TRAVELED 


b) MIDCHORD PITCHING MOMENT COEFFICIENT 

Fig. 3 Effect of far-field boundary dislancc on and for step 
change in angle of allack. 


Table 1 Indictal response to a step change in angle of attack: 
NACA 64A006, a = 1 deg and are given In units of chord) 

Boundary condition 



Reflecting 

Nonreflecting 
(1st approx) 

A/. =0.80 

Mesh 

103x81 

93x71 

J'ma, 

25 X 79 c 

3.8x21. 5c 

CPU' 

1.19 

1.00 

II 

p 

bo 

Mesh 

103x81 

93x79 


25 X 79 c 

3.8x6lc 

CPU* 

1. 10 

1. 00 


"CPU lime is normalized by the nonreflecting case. 


Table 2 Sinusoidally oscillating trailing-edge flap: 
NACA64A006,tf = slnw/ 



Boundary condition 



Nonreflecting 


Reflecting 

(1st approx) 

Mesh 

93x77 

93x61 


3.8x47c 

3.8x7.6c 

CPU* 

1.24 

1.00 


"CPU lime is normalized by the nonreflecting case. 







1406 


D. KWAK 


AIAA JOURNAL 


TIME, CHORD LENGTHS TRAVELED 



I '0.18 


o 



t = 0.18 


£> 


£... '..i 

t = 0.18 



a) LARGE GRID SOLUTION 


b) REFLECTING B.C. 


Fig. 6 Pressure contour for an impulsively started NACA 64A006 airfoil from rest. 


cj nonreflecting B.C. 


LARGE GRID SOLUTION 



TIME. CHORD LENGTHS TRAVELED 

Fig. 7 Upper surface pressure coefficient for an impulsively slarted 
NACA 64A006 airfoil. 

As shown in Fig. 5, when the reflecting boundary is located 
loo close to the airfoil, reflected waves cause the shocks to 
move slightly downstream in both the type A and type B 
motion. The reflecting boundary conditions on the top and 
bottom, =0, model solid-wall boundaries. Hence, this 
close-in imposition restricts the flow and produces stronger 
shocks. However, applying nonreflecting conditions, waves 
incident on the far-field boundaries are mostly absorbed, 
yielding reasonably good results as shown in Figs. 4 and 5. It 
is to be noted that slight changes in profile (Fig. 5) produce 
significant changes in unsteady loads (Fig. 4). The gain in 
computational efficiency, due to the first approximate 
conditions as well as the reduction in computational domain, 
is shown in Table 2. 


C. Impulsively Started Airfoil 

To visualize the behavior of waves originating from the 
airfoil and then propagating throughout the computational 
domain, an NACA 64A(X)6 airfoil is impulsively slarted from 
rest at time zero with A/* =0.85. This is equivalently done in 
LTRAN 2 by turning on the freestream from rest. To see the 
reflection more clearly, the far-field boundaries are placed 
very close to the airfoil = 1.22 c, = 1.08 c). The 

propagation of disturbances is demonstrated by the pressure 
contour of the upper half plane of the computational domain. 
In Fig. 6, the results are shown at three sequential limes 
measured by chord lengths of airfoil traveled. 

Figure 6a is obtained from the large grid solution by 
looking through the window of 1 .22 x 1 .08 c. In this sequence 
taken from a movie produced from the calculation, waves 
from the airfoil propagate outwardly only during the time 
observed as though the domain were infinite. In Fig. 6b, the 
old reflecting conditions are used. Physically, this is 
equivalent to placing solid walls on the boundaries. As time 
increases, the influence of the boundaries, especially the y 
boundary, becomes more apparent and, at / = 7.1, the flow is 
shown to be choked. This illustrates a possible influence that a 
solid wind-tunnel wall can have on experimental results. 
When the nonreflecting conditions arc used (Fig. 6c), the 
pressure contour map more closely resembles the large grid 
solution. As shown by the upper surface pressure coefficient 
plot in Fig. 7, the reflecting conditions contaminate the near- 
field solution much faster. Since the nonreflecting conditions 
used are simple approximations to perfectly absorbing 
conditions, some reflections still exist. 

V. Conclusion 

The first approximate nonreflecting far-field boundary 
condition procedure of Engquisl and Majda has been suc- 
cessfully applied in solving various practical problems using 
the unsteady small-disturbance transonic flow code, 
LTRAN2. (As discussed in Sec. IIB, their first approximation 
is a form of the characteristic equation. Therefore, the 
boundary condition procedure based on the first ap- 
proximation alone may be regarded as a characteristic con- 
dition procedure.) It has been demonstrated successfully that 
the computational domain can be reduced considerably using 
their boundary condition procedure. However, since the mesh 
system is stretched in LTRAN2, the number of computational 
grid points is not linearly proportional to the distance of the 
boundary from the airfoil. The gain in computer time after 
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implementing the new conditions is problem dependent and, 
for the cases tested using the first approximate conditions, the 
time saved is 10-24*^7o. Further reduction in the computational 
mesh may be achieved by higher-order approximations. 
However, the benefit obtained by the reduced number of 
mesh is overshadowed by the increased computing time 
required to implement these high-order boundary conditions. 
Since the approximate conditions applied here are very simple 
and do not permit any waves to propagate in from the out- 
side, modifications can be made to accommodate the waves 
which should come back into the domain. However, 
especially in the production version of a code, the simplicity 
of the first approximate nonreflecting conditions and the 
computational efficiency gained are more than adequate for 
applications in the design process. 
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